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, Abstract. We extend to several dimensions the result of K. Seip 

and Y.I. Lyubarskii that proves the existence of Riesz basis of 
i exponentials for a finite union of intervals with equals lengths. 



^ . 1. Introduction 

' Let E G M.'^ he a bounded set with no empty interior. We consider 

■ the Paley- Wiener space of functions band-hmited to E 

PWe = {F e L^iW^) ■ supJ^F C E}, 

endowed with the L'^{W'') norm. Here we denote as 



(N 



the Fourier transform of F. 



E 



O . By the Paley- Wiener theorem PWe is a reproducing kernel Hilbert 

Q ! space of entire functions with kernel 

1: 

S ■ Then, for any F G PWe, 

> 

X 
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We say that a sequence A = {A„} C M'^ is sampling for PWe if 
there exist constants A,B>0 such that 

(1) a[ \F{x)\^dx<y^\F{K)\'<B f \F{x)\'dx, 

for all F e PWe. 

We say that a sequence A = {A„} C M.'^ is interpolating for PWe 
if for each {a„} G £^ there exists a function F G PVTe such that 

(2) F(A„) = a„, n G Z. 

If A is both sampling and interpolating we say that it is a complete 
interpolating sequence. 

Sampling and interpolation sequences in Paley- Wiener spaces pro- 
vide a mathematical model of stable recovery and data transmission 
in signal theory regardless of the computation method. For the mean- 
ing of this concepts we refer to jLan67bj . In 1967 H. J. Landau found 
necessary conditions for the existence of such a sequences in terms of 
some densities, as a consequence a complete interpolating sequence 
must have density equal to the Landau-Nyquist rate, see Theorem 

The aim of this paper is to find a complete interpolating sequence 
for PWe, where E is a. finite union of mutually disjoint half-open 
cubes 

p d 

(3) ^ = U C R'', where Qj = a'j + p), /3 > 0. 

j=l s=l 

Our result allows us to construct, for any bounded E with no empty 
interior, sampling or interpolating sequences with density arbitrarily 
close to the critical one, see Corollarv l3.4[ This implies that Landau's 
results |Lan67j about densities can't be improved. 

Moreover, we find a recovery process in PWe whose stability is 
assured by the sampling property, see Theorem 13.31 This recov- 
ery process can be considered as a shifted version of the Shannon- 
Whittacker-Kotelnikov theorem. 

The method used is a generalization to Mf^ with c? > 1 to that of 
|LS97j for M. It is worth noting that the reconstruction formula (fT^ . 
which comes from A. Kohlenberg's work |Koh35j . has been used by A. 
Faridani to produce sampling schemes for computerized tomography 
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|Far9nj . Our result provides the stability that their result lacks. In 
fact, Faridani proved a more general version of the reconstruction 
formula ()14j) that holds for smooth functions with spectrum in any 
compact set of R'', and there is a version even for the group context 
|BF02j . But in such a general context there is no hope of stability. 

Finally, we observe that we can define for 1 < g < oo the Paley- 
Wiener space PW^ of entire functions given by 

n^) = j^2lj(iy'-'dt, zee' 

with f E L"^ (E), for q = q/{q~ 1). With the natural modifications 
the main result (Theorem I3.3|l holds also in PW^ for 1 < g < oo 
because the projection of L'^iW'-) into PWi^E) is bounded and we 
have also a Plancherel-Polya type inequality j Nik75t 3.2.2.]. 



2. Preliminaries 

In this section we review some results and definitions that will be 
used later on. 

Let V = {vn} be a sequence in a separable Hilbert space H. We say 
that f is a frame for H if there exist constants A,B > such that 

A\\u\\^ <^\{u,Vn)\^ < BWuW^, ueH. 

n 

We say that v = {f„} is a Riesz sequence or, equivalently, that it is a 
Riesz basis in its closed linear span, if there exist constants A,B > 
such that for every finite scalar sequence {c„} one has 

Now we can translate our problem concerning sequences into the 
basis language. A sequence A = {A„} C Mf^ is interpolating for PWe 
if and only if the sequence {K{-, A„) / a/ K{Xn, A„)} is a Riesz sequence 
for PWe, and this is also equivalent to the assumption that 

^(A) = {e'^--^ : G A} 

is a Riesz sequence for L'^{E). Similarly, A = {A„} C M'^ is sampling 
if and only if {K{-, A„,)/ a/ K{Xn, A„)} is a frame for PWe, and if and 
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only if 

^(A) = {e^^™-^ : A„ G A} 

is a frame for L'^{E). Then A is both samphng and interpolating if 
and only if S{A) is a Riesz basis for L'^{E), i.e. S{A) is the image of 
an orthonormal basis by a bounded invertible operator, see |LS97j . 
We say that a sequence A = {A„} C is uniformly discrete if 

inf \Xj — \k\ > 0. 

This infimum is called the separation constant of A. 

For h > and x = {xi, . . . , Xd) denote by Qh{x) the cube: 

d 

Qh{x) = Y[[xs-h/2,Xs + h/2). 

s=l 

Then {Qh{hm)}^^'id is a disjoint cover of M'^ for any h > 0. 

Let A = {A„} C M"' be a sequence. Following Beurling we define 
the upper uniform density D~^{A) and lower uniform density D~{A) 
of A as follows: 

D+(A) = hmsup^^^ and ^-(A) = liminf 

h^oo 1^^°° h"' 

where 

n+{h) = sup #(A n (5/,(a;)) and n'{h)= inf #(A n (5/,(a;)). 



If this densities coincide we say that A has uniform density and we 
write D{A) = D+{A) = D-{A). 

The following result, which provides necessary conditions for sam- 
pling and interpolation, is due to H. J. Landau |Lan67j . 

Theorem 2.1. If A is an interpolating sequence for PWe then A is 
uniformly discrete and 

D+{A)< 



If A is a sampling sequence for PWe there exists A' C A uniformly 
discrete and sampling subsequence such that 

D-(A') > 
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Then Riesz basis can occur only when D{A) = |£'|/(27r)'^. This 
particular value is called the Landau-Nyquist rate. 

It is quite easy to show that no characterization can be obtained 
for sampling and interpolation in several dimensions, nor in the one- 
dimensional multi-band case, using only these densities. 

Unlike the case of a single interval, where there is a complete char- 
acterization for the complete interpolating sequences |Pa79j . little 
is known for the same problem with several intervals or in several 
dimensions. Let us comment some of the known results. The reason- 
able question in this more general setting is to find some complete 
interpolation sequence rather than to obtain a complete characteri- 
zation. For a finite union of commensurable intervals, for two general 
intervals or for particular configurations with more than two inter- 
vals, it is known that such sequences exist |LS97j . In there are 
also complete interpolating sequences for PWe if -E is a convex poly- 
gon symmetric with respect to the origin |LR00j . In M°' for d > 2 



the only examples, as far as we know, concern the spaces PWe with 
parallelogram domains E. For more on this problem see |Sei04j . 

Our result provides the first example of non-convex and non-symmetric 
set with Riesz basis of exponentials in several dimensions. This re- 
sult contrast with the more rigid analogue problem of orthonormal 
basis of exponentials for L'^{E), where it is known that symmetry is 
necessary |KolOOj . 

Our result allows us to construct, for any E G M."^ bounded with 
no empty interior, sampling or interpolating sequences for PWe with 
density arbitrarily close to the Landau-Nyquist rate, see Corollary 
13.41 Unfortunately this method doesn't seem to be useful to find 
complete interpolating sequences for PWe with E arbitrary. 

Next we will show that the separation property of a sequence is 
equivalent to the right inequality in (Q). This is a well known result 
but we include its proof for the sake of completeness |You80l pp 82- 
83]. 



Lemma 2.2. Let E G be a bounded set with no empty interior 
and A = {A„} G be a sequence. A is a finite union of uniformly 
discrete sequences if and only if there exists a constant B > such 
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that 

(4) 5^|F(A„)|'<5 / 

for all F G PWe- Moreover, the constant B > depends only on the 
diameter of E and on the separation constant of A. 



Proof. Suppose, to reach a contradiction, that we have (jH) but A is 
not a finite union of uniformly discrete sequences. We can suppose 
that G int-E and then there exists an e > such that <5e(0) C E 
and 

sup #(A n Qhihm)) = oo, 

where h = ^. Choose a sequence {mjjjgz C Z'^ such that 
#(A n Qhihruj)) oo, j ^ oo. 

The functions fj{x) = e~^^"^^'^XQe{o)i^) have support on E and 

1 



d 



fj{x)e"'''dx 



are in PWe- Now 



inf inf > 0, 



because 



inf \Fj{z) 



> 



-, — 73777 inf 

(27r)'^/2 .6Q,(0) 

e 



n 

s=l 



smez 



inf 

277 |t|<l/2 



sin 77t 



77t 



and therefore 



l^.|A| 



£2 ^ oo, J ^ OO, 



but IIF,- 



1/2 
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To prove the converse, we can suppose that A is uniformly discrete, 
i.e. there exists r > such that 

|Aj - Afcl > r > 0, j ^ k. 

Take an interval that contains E 

d 

E C JJ[-a^a'^]. 

s=l 

Then F G PWe is of exponential type v = {u^, . . . , u'^), where < 
u'^ < , s = 1, . . . ,p. From the mean- value property one has, for all 
Ps > 0, 

^(^) = 7^ [ ^(^' + Pie'"', + pde''^)de, X G R'^. 

Integrating on [0, 5]°', for some 5 > 0, 

= 7^ / F{z)dm{z), 
\2 J (27r)^ Jj^d^^ 5(^.5) 

where i?(a;'',5) = {z & C : {x"^ — z\ < 6}. Then by Cauchy-Schwarz 
inequality 

\F{x)\'<b[ if \F{i + ir,)\^dO\dr^, 

where i? > depends on 5 and d. 

For all F entire function of exponential type (z/^, . . . v'^) such that 
F G L^(M'^), we have the following inequality of Plancherel-Polya 
type 

(5) ||F(- + ^r7)|U.(i,.) < e^'-i'^=l''^l||F|U.(K.), 

see |Nik75[ 3.2.2.]. 
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Finally, if we take 5 = ^ the cubes Q2s{^j) sue pairwise disjoint. 
This together with (0), gives 

5^|F(A,)|^ < B [ Y.I \F{^ + ^r^)Mdr^ 

< B [ ( [ \Fi^ + tr])\'dAdr^ 

■ 

Following Kohlenberg, |Koh35j . we define the sampling procedure 
of order p for a function F G PWe such that J-'F G C^{E), as follows 

p 

(6) Git) =J2Y1 ^(^^ + ^')'^'(* -Wn- h). 

1=1 n&d 

Here ly is a real non-singular d^ d matrix, ki G M*^, and Si G PWe- 
Since {Vrn + A;i}„gzd is uniformly discrete (when the ki are different) 
and F tends to zero faster than any polynomial this series converges. 
We want to determine ki and Si in order to get G = F for all 
F G PWe. 

First we compute the Fourier transform of this sampling procedure. 

Proposition 2.3. Let F, Si G PWe, f = J^F e C^{E), si = J^Si 

and G given by Then 

(7) 

' ' 1=1 raeZd 

Proof. This result follows easily from Poisson's Formula and it can 
be found in |Far90j . Consider the Fourier transform of G 

p 

J^G{u) = J^sKcc') J2 F{Wn + ki)e-'''<'^''+^'\ 

1=1 neZ<* 
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Poisson's formula gives 

where F^, = F{- + ki). m 

3. Construction and main result 
Let E he a. finite union of lialf-open cubes 

p d 

E=[jQjC where Qj = a'j+(3), (3 > 0. 

j=l s=l 

We want to determine the functions and shifts Si{- — ki), I = 1, . . . ,p 
in order that the reconstruction formula 

p 

F{t) = ^(^^ + ^')'^'(^ -Wn- ki) 

holds for all F G PWe such that f = J^F e C^{E). We know, by 
Proposition 12.31 that for such functions this is equivalent to 



=1 



From now on we take W = ^Id and, as / = off of £" we define 
si{uj) = if C(j ^ i?. To determine si on E we split E in the parts 
which are overlapped by one or more of the sets 

E + 2Tr{W-yn, n G Z"^. 

Let us see how this partition works. 

Consider two cubes Qj and Qk- There exists a unique rijk G Z*^ 
such that 

ijk = c^k- Pfijk eQj, 
where ak = (a^, . . . , af) is the lower left point of Qk- Observe that 
Ujj = 0. Fixed j G {1, . . . ,p} split Qj into the half-open intervals 
formed by intersection of Qj with the 2'^ quadrants around 7jfc, and 
do the same for all of = 1, . . . , p. Thus, we split each cube Qj in this 
way in at most p'^ intervals (this can be easily proved by induction, 
projecting a face of the cube into M^~^) which we denote QjS = 
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Figure 1. 

1, . . . jP'^, the subscript indicates the cube in which the rectangle is 
contained. From now on we assume that we have p'^ intervals on 
each cube Qj. The minor changes needed if that is not the case are 
left to the reader. We have for each j = 1, . . . ,p 

pd 
s=l 

Although 7jfc + [3njk = ak E Qk and jjk G Qj, it can happen that 
for some Qj C Qj 

iQ'^ + (3n,k) n Qfe = 0, 
see figure ^ To solve this difficulty we choose a correction 

C,,(A;)G {-1,0,1}^ 

in such a way that 

Q- + f3{n,k + C,s{k)) = Q'k, 

for some G {1, . . . ,p'^}. 

This shows that the intervals that form the cubes are always the 
same, but put in different orders, i.e. for any Ql C Qi we have a 
unique rectangle in Qj that is translated by jSZ'^ in Qf. We call this 
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rectangle Qj- Then the intervals on each Qj are named depending on 
their translate on Qi- Using this convention we have for s = 1, . . . jp'^ 

Aj — X 'p 

Q{ + f3{n^k + Cu{k))=Ql, 

and 

(9) Ukj + Cks{j) = nij - Tiik + Cis{j) - Cis{k). 

This is so because the translation by f3Z'^ that moves Qj to their 
counterpart in Qk is unique. 
Let 

Cjs = {njk + Cjs{k)}k=i,...p 
and notice that = njj + Cjs{j) G Cjs and as the cubes are disjoint 
the elements in Cjs are all different. 
Now we have, for u ^ Qj 

uj + (3n e E iff n e Cjs, 

then for a; G equation (jSJ becomes 

d P 



J2 fiu; + (3n)f2^ii^y^'^' 



(3 



(10) 



271"/ , , 

If we find for all j and s a solution xi = xi{j, s) I = 1, . . . ,p of the 
linear system 

ELia^^e^'""''' =0, if meC,s\{0} 
Er=i = , otherwise 

we have finished, because defining si{lj) = xi{j, s) for G and 

Si(^) = EE^'(^'' 

we have (jH)). 

We claim that there exist ki G M'^ for / = 1, . . . ,p such that (fTUj) 
has a solution for all j = 1, . . . , p and s = 1, . . . ,p'^. 



12 JORDI MARZO 

Indeed, letting zis = e^^^i for / = 1, . . . ,p, s = 1, . . . ,d and zi = 
{zii, . . . , zid) the linear system in (fTUI) has solution if and only if there 
exists zi G T'^, I = 1, . . . ,p such that 

j=l s=l 

Indeed, letting 

mjs = ( niin nj^ + C],(fc), . . . , niin n'^j,^ + Cf^ik)), 

k=l,...,p ■' k=\,...,p 

then 

^^^[.X^l Sk,l=l,...,p) — rrijs rrijs 

where poljs is a non-zero polynomial in C^'^. Now the first column 
expansion of the determinant gives us 



k=l 



Since the elements in Cjs are different, an induction reasoning pro- 
vides the conclusion. 

Moreover, the zero set of V has zero Lebesgue measure in the 
polydisc JP'^, because T> is an entire function |R,on74| p. 218]. Then 
for almost all {ki, . . . ,kp) G [0, ^)^'^, and all j and s the linear system 
in (fTUI) has a solution. 

So far we have shown that for F G PWe such that J^F G C^{E), 
I 



there exist ki G M.'^ with 



(12) ©(e'^'^V . . , e*^'^^') ^ 0, 

so that taking 



p p" 



(13) Slit) = :F-\si{t) =Y,Y.^i{3.s):F-\Qs{t), I = 1 

j=l s=l 

we have 

p 

(14) Fit) = Y^^F{Wn + ki)Si{t-Wn~ki). 

1=1 neZt* 
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The following results allows us to extend this formula to all F E 
PWe- 

Lemma 3.1. Take ki G M"' such that il^) holds and functions Si as 
in 173]). Then, if a = Uf^^a', a' G i"^, the series 

p 

(15) F[a]{t) = J2Y1 "^nSiit -Wn- k), 

1=1 n&Zd 

converge in L'^{M.'^). Moreover F[a] is an infective operator: there 
exists B > independent of a such that 

(16) \\F[a]\\L2(^Rd) < B\\a\\p. 

Corollary 3.2. Taking functions Si as in il!^) . and ki G M'^ such 
that / lijjj) holds, then 

Si{Wn + ks- k) = 5nfi5i^s, l,s = l,...,p. 

The proofs are given in the next section. We can now prove our 
main result. 

Theorem 3.3. Let E dW^ he a finite union of cubes as in (0). // 
W = j-Id and h G W are such that (C^ holds then 

K{k^,...,kp) = U^i^^{Wn + ki}^^^, 

is a complete interpolating sequence for PWe, i-e. there exists con- 
stants A,B>0 such that 

p 

(17) A||Ff <^ J]|F(iyn + A;0|'<5||Ff, F e PWe. 

1=1 neZd 

Given any a = Uf^j^a' with a' G i"^ there exists a unique F G PWe 
such that 

(18) F{Wn + ki)=al, neZ'^,l = l,...,p. 

Moreover, taking functions Si as in il^) . the reconstruction formula 
l{14\ ) holds for all F G PWe, with convergence in L'^—norm and uni- 
form convergence on products of horizontal strips. 
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Proof. By the preceding construction we know that (fT^ holds for all 
F G PWe such that J-'F G C^{E). Now by Lemma EH] the operator 
on PWe defined by 

p 

F I — ^ XI 5Z ^(^^ + -Wn- k), 

is bounded and the sampling formula ()14|) holds for PWe- 

Since A(A;i, . . . , fcp) is uniformly discrete we have the right inequal- 
ity in ()17p . In order to prove the left inequality, notice that by Lemma 
13. II we have that '^n^ii' — Wn — ki) converges for all a' G i"^, there- 
fore there exist a constant B > such that for F G PWe 



Ln 



\{F,Si{--Wn-ki))\' <B\\F\\', 
see |Chr03| Lemma 3.2.1.]. Then for F G PWe we have 

Si{--Wn-ki))\' 



\ l,n / \ l,n 

< B\\Ff (^\F{Wn + k)\^^ 



Defining for a = Uf^^^a', a' G F[q;] G PWe as in (fT3)) . and using 
Corollarv 13.21 we get the interpolation property (fTH|) 

p 

F[a](iym + A;,) = ^ ^ a^5;(iy(m - + A;, - A;;) 

p 

Finally, the uniform convergence on products of horizontal strips 

{zeC':\^z^<C,,j = l,...,p}, 
follows from the Phragmen-Lindelof inequality, see |You80| p 84] . ■ 



RIESZ BASIS OF EXPONENTIALS FOR A UNION OF CUBES IN R'' 15 



Corollary 3.4. Let E dW^ he a hounded set with no empty interior. 
For all e > there exist uniformly separated sequences A^, C M.'^ 
such that Ae is sampling for PWe, A^ is interpolating for PWe and 

< {2tiYD{K,) - 1^1 < e, < 1^1 - {2-kYD{K') < e. 

Proof. It is enough to choose E^, E'^ as in Q such that E'^ G E G E^ 
and \Ei,\ — \E\, \E\ — \E'^\ < e. Now we take a complete interpolating 
sequence A^ for PWe^ (that obviously is also interpolating for PWe) 
and Ae a complete interpolating sequence for PWe^ (that is also 
sampling for PWe) and apply Landau's Theorem 12.11 ■ 



In this section we will proof the technical results used in proving 
Theorem 12.31 

Proof. [Lemma 13. Ij Consider, given M e N 



4. Proofs of technical results 



p 



FM[a]{t) = J2Y1 ""nSiit -Wn- h). 



1=1 \n\<M 



For all / G 



{W^) we have 



(FuHf) = J2J2x,ij,s)J2aUe-<'^-^'^^ 



XQ^f) 



l,j=l s=l \n\<M 



p p"^ 





l,j=l s=l \n\<M 
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where Vpwosif) = ^^^{Xq^^ f) is the orthogonal projection of / G 
j 

LF'{W^) into PWqs. Then, using Lemma f2.2l {<£). we obtain 

\{FM[a]J)\ < ^)|||a'||,.rp^,.(/)IU^ 

l,j=l s=l 

p p p'^ 

i=l j=l s=l ^ 

p 

< i?||/|U2^||a'||,2, 

where the constant B does not depend on a. 

Now we want to prove the injectivity part. By taking Fourier 
transform in ()15|) we obtain 



j=l s=l 1=1 KneZ'i 

We call 

If we show that 

fi{uj) = 0, l = l,...p, 

we are done. The functions // are /5— periodic on each variable, hence 
it is enough to proof fi\Qi = 0. We will prove this on each Ql sepa- 
rately. The key is property Q. 

Let us take u E Q^. Then, as fi{u + [3{nik + Cis{k))) = fi{uj) 
and uj + /?(nifc + Cisik)) G Qks, substituting into (|T9|l we obtain, for 
A; = 1, . . . ,p 

p 

1=1 
p 
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We have a linear system with unknowns e^'^^'^'- fi{uj) and coefficient's 
matrix 

We will see that this system has only the trivial solution by showing 
that it has invertible coefficient's matrix. This point becomes quite 
easy with the notation we have adopted. Indeed, the linear system 
(fTUI) that we have used to determine si in Q\ is given by 

Now by the elements of the matrix product are 
p p 

1=1 1=1 
p 

1=1 

which is the desired conclusion. H 

Proof. [Corollary IH.2j Let s G {l,...,p} be fixed. The sampling 
result applied to Ss{- — h) E PWe yields 

p 

Ss{t - A;) = 5^ 5^ ^,(H^n + A;; - h)Siit -Wn- h). 

1=1 neZ'' 

Now applying Lemma f3. II to 

p 

1=1 neZd 

we obtain the result. ■ 




References 

[Far90] A. Faridani, An application of a multidimensional theorem to computed 
tomography, Contemp. Math. 113, 65-80, 1990 

[BF02] H. Behmard, A. Faridani, Sampling of handlimited functions on unions 
of shifted lattices, J. Fourier Anal. Appl. , vol. 8, 1, 2002 

[Koh35] A. Kohlenberg, Exact interpolation of Band-Limited Functions, 
J.Appl.Phys., 24, 12, 1432-1436, 1935 

[KolOO] M. Kolountzakis, N on- symmetric convex domains have no basis of expo- 
nentials, Illinois J. Math., 44, no. 3 , 542-550, 2000 



18 



JORDI MARZO 



[Lan67] H. J. Landau, Necessary Density Conditions for Sampling and Interpo- 
lation of Certain Entire Functions, Acta. Math., 117, 37 52, 1967 

[Lan67b] H. J. Landau, Sampling, Data transmission, and the Nyquist rate, Proc. 
IEEE, voL 55, no. 10, 1701-1706, 1967 

[LROO] Y. I. Lyubarskii, A. Rashkovskii, Complete interpolating sequences for 
Fourier transforms supported by convex symmetric polygons, Ark. Mat. 
38, 139-170, 2000 

[LS97] Y. I. Lyubarskii, K. Seip, Sampling and Intepolating Sequences for 

Multihand- Limited Functions end Exponential Bases on Disconnected Sets, 

J. Fourier Anal. AppL, 3, 5, 597-615, 1997 
[Nik75] S. M. Nikol'skii, Approximation of Functions of Several Variables and 

Imbedding Theorems, Springer- Verlag, Die Grundlehren der mathematis- 

che Wissenschaften in Einzeldarstellungen, 205, 1975 
[Chr03] O. Christensen, An introduction to frames and Riesz basis, Birkhauser, 

ANHA, 2003 

[Pa79] B. S. Pavlov, The basis property of a system of exponentials and the con- 
dition of Muckenhoupt, Dokl. Akad. Nauk SSSR, 247, no. 1, 37-40, 1979 

[Ron74] L. L Ronkin, Introduction to the theory of entire functions of several 
variables, Trans. Math. Monographs, AMS, Providence, RI, 1974 

[Sei04] K. Scip, Interpolation and sampling in spa,ces of analytic functions. Uni- 
versity Lecture Series, 33. American Mathematical Society, Providence, 
RI, 2004 

[You80] R. Young, An introduction to nonharmonic Fourier series. Academic 
Press, New York, 1980 (revised first edition 2001) 

Departament de matematica aplicada I ANALisi, Universitat de Barcelona, 
Gran via 585, 08071 Barcelona, Spain 
E-mail address: jmarzo@mat.ub.es 



